Theory of magnetic switching of ferroelectricity in spiral magnets 
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We propose a microscopic theory for magnetic switching of electric polarization (P) in the spin- 
spiral multiferroics by taking TbMnOa and DyMnOa as examples. We reproduce their phase dia- 
grams under a magnetic field J?ox by Monte-Carlo simulation of an accurate spin model and reveal 
that competition among the Dzyaloshinskii-Moriya interaction, spin anisotropy, and spin exchange 
is controlled by the applied //ox, resulting in magnetic transitions accompanied by reorientation or 
vanishing of P. We also discuss the relevance of the proposed mechanisms to many other multifer- 
roics such as LiCu202, MnW04, and Ni3V204. 

PACS numbers: 77.80.Fm, 75.80.+q, 75.30.Gw, 75.47.Lx 



Concurrently magnetic and ferroelectric materials, i.e. 
multiferroics, offer prospective systems to attain mag- 
netic control of electricity via magnetoelectric (ME) cou- 
pling [HQ- It was experimentally demonstrated that 
an external magnetic field (i?ox) can cause reorientation, 
emergence, and vanishing of ferroelectric polarization P 
in many spin-spiral multiferroics such as /JMnOa (^=Tb, 
Dy, Eui_^Y^, etc) [l-IHl, LiCuaOz 1, MnW04 B: and 
Ni3V204 These ME phenomena arc currently at- 
tracting enormous interest, and a thorough understand- 
ing of their mechanisms is an urgent issue. However, the 
number of theoretical studies is very few despite many 
experimental reports. Naively, the applied ii/ex can de- 
termine the direction of P by controlling the conical spin 
structure via Zeeman coupling, but there are many ex- 
amples that do not obey this simple picture. 

In the spin-spiral multiferroics, inherent spin frustra- 
tion as an origin of the spiral magnetism inevitably re- 
duces the spin-exchange energy, and hence increases the 
relative importance of other tiny interactions, e.g. the 
single-ion spin anisotropy and the Dzyaloshinskii-Moriya 
(DM) interaction. Consequently, the magnetic switch- 
ing of P in this new class of multiferroics is governed 
by their fine energy balance tuned by iifcx, which cannot 
be understood from a simple interplay between Zeeman 
coupling and the spin exchanges. 

In this Letter, by taking the Mn perovskites TbMnOa 
and DyMnOa as examples, we propose a microscopic the- 
ory for the magnetic control of P in the spin-spiral multi- 
ferroics. Their puzzling T-H^x phase diagrams are repro- 
duced by the Monte-Carlo (MC) analysis of an accurate 
spin model. Our microscopic theory reveals that the ap- 
plied iifex controls conflicts among the spin exchanges, 
spin anisotropy, and DM interaction, resulting in mag- 
netic transitions accompanied by reorientation or vanish- 
ing of P. The mechanisms proposed here are relevant to 
many other spin-spiral multiferroics such as LiCu2 02 @, 
MnW04 0, and Ni3V204 0. We also discuss the in- 
fluence of effective magnetic flelds from rare-earth / mo- 
ments. 



The ferroelectricity in these materials is described by 
the spin-current model 0, as given by P oc Q x x, 
where Q is a propagation vector of the spiral and x oc 
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Si X Sj is the vector spin chirality. As shown 
in Fig [1Kb), the Mn spins in TbMnOs and DyMnOs are 
rotating within the be plane (x||a) to form a transverse 
spiral with Q\\b [ll|, and thus P\\c is realized. 

In Figs. Hc)-(f), we briefly summarize the puzzles in 
.RMnOa [4|. The applied H^x induces the magnetiza- 
tion M\\Hex via Zeeman coupling, and hence forces the 
spin structure to be conical where xll-f^ox- When we 
apply HcxWQ [see Fig. IHc)], we expect a longitudinal 
conical spin order with xHQ. In this case, P should 
be zero within the spin-current model. Thus we ex- 
pect vanishing of P when we apply Hcx\\b {Phnm set- 
ting) to TbMnOs and DyMnOa. However, reorienta- 
tion of P from P\\c to P\\a is observed in reality [see 
Fig. [He)]. A neutron-scattering experiment conflrmed 
that this P reorientation results from the spin-chirality 
flop from x||a to x||c [l^, [l^. This discrepancy has 
been naively attributed to the influence of / moments 
on the rare-earth ions thus far 15, loj. However, a sim- 



ilar behavior has been observed also in LiCu202 with- 
out / moments @, suggesting an essentially new mech- 
anism. Mostovoy reproduced the flop by introducing 
higher-order anisotropics in a phenomenological theory 
although their microscopic origins are unclear On 
the other hand, the application of Hcx-^Q is expected to 
stabilize a transverse conical spin order with X-LQ. As 
shown in Fig.[IJd), we expect the a&-plane transverse con- 
ical order with P\\a when we apply Hcx\\c to TbMnOa 
and DyMnOa. However, in TbMnOa, the first-order 
transition to paraclcctric {P=0) phase is observed under 
Hcx\\c as shown in Fig.[TKf). The ifex-induced vanishing 
of P is also observed in MnW04 and Ni3V204 01 . 

To solve these puzzles, we start with a classical Heisen- 
berg model on a cubic lattice, in which the Mn 5=2 
spins are treated as classical vectors. The Hamiltonian 
is given by 7^ = Tij + Hsia + Hum + "Hzooman- The first 
term Hj = E<i j> Jij'^i ' '^j describes spin-exchange in- 
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FIG. 1: (color online), (a) Crystal structure, spin exchanges, 
and local axes ^i, rji, and in RMnO^. Here FM (AFM) de- 
notes (anti)ferromagnetic exchange, (b) fee-plane transverse 
spin spiral in TbMnOs and DyMnOa, which induces ferroelec- 
tric polarization P\\c. (c) [(d)] Application of iifcx||b [i?cx||c] 
is expected to stabilize the longitudinal [transverse] spin spi- 
ral with magnetization M\\Hcx where P=0 [P\\a] is expected 
within the spin-current model. (e)[(f)] Experimental T-Hc^ 
phase diagram of DyMnOa [TbMnOs] for Hc^\\b [Hc^\\c] from 
Ref. 0], which shows reorientation of P from P||c to P||a 
[disappearance of P] . Here FE (PE) denotes ferroelectric 
(paraelectric) phase. 



teractions as shown in Fig. [TJa). The second term Hsia, 
denotes the single-ion spin anisotropy, which consists of 
two parts as H^ia = n^.^+n^^ with n^^^ = S^, and 

= EY.^{-lr''^'^{Sl - S^). Here 6, m and C, are 
the tilted local axes attached to the ith MnOe octahe- 
dron . The term T-L^g^ causes the hard-axis anisotropy 
along c, or, equivalently, the easy-plane anisotropy in the 
ab plane. The third term "Hdm = J2<i j> ' i'^i ^ ■^i) 
represents the DM interaction where the vectors dij are 
defined on the Mn(i)-0-Mn(j) bonds, and are expressed 
by five DM parameters, aab, Pab, lab, etc, and /3c [I!]. 
This term consists of two parts, W^yi and T-L^yi-, where 
^DM (^dm) associated with the DM vectors on the 
in-plane (out-of-plane) Mn-O-Mn bonds. The last term, 

/ccman = .g/iB Yji Si ■ -f^in; stands for the Zeeman cou- 
pling. Note that the Mn spins feel the internal magnetic 
field i?in, which consists of two contributions, i.e., the 
applied field i?ox and the effective field Hfd from the 
/ moments. This model has successfully reproduced the 
phase diagrams of i^MnOs at Hc^=0 [l9|. 

We have microscopically determined the values of Jab, 
Jb, Jc, D, and E, and have estimated the values of five 
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FIG. 2: (color online). Theoretical T-Hin phase diagrams of 
(a) TbMnOs and (b) DyMnOs for Hin\\b. (c) T profiles of 
specific heat C{T) and spin chiralities x-yi^) il—'^, ^, c) for 
TbMnOs at Hl^=& T. (d) Spin structure in the &c-plane spiral 
state at i/in=0, and arrangement of the a-axis components of 
DM vectors on the out-of-plane Mn-O-Mn bonds. The sym- 
bols and ® express their signs, i.e., positive and negative, 
respectively. In the inset, the arrows (dashed lines) show the 
spin directions in the presence (absence) of DM interaction. 



DM parameters in Ref. [l^. We perform calculations 
using two sets of the model parameters (A and B) as 
(A) {Jab, Jb, Jc)=(-0.74, 0.64, 1.0), {D, E)={0.2, 0.25), 
(aab, Pab, 7afc)=(0.1, 0.1, 0.14) and (a^, /3c)=(0.48, 0.1), 
and (B) {Jab, Jb, Jc)=(-0.7, 0.99, 1.0), {D, S)=(0.22, 
0.25), {aab, Pab, 7afc)=(0.1, 0.1, 0.14) and (a^, /3c)=(0.45, 
0.1). Here the energy unit is meV. These parameter sets 
give the 6c-plane spin spirals propagating along the b axis 
with wave numbers (5b=0.37r and Qb=0.47r, respectively. 
They reproduce well the spiral spin states in TbMnOa 
(Q6=0.287r) [nj and DyMnOg (g6=0.397r) [| at i7ox=0. 

We analyze this model using the replica-exchange MC 
technique [20[ . Each exchange sampling is taken after 400 
standard MC steps. After 600 exchanges for thermaliza- 
tion, we typically perform 1000 exchanges for systems of 
7V=40x40x6 sites with periodic boundaries. 

In Figs, mja) and (b) we display theoretically ob- 
tained T-Hia phase diagrams of TbMnOa and DyMnOa 
for i?i„||6, respectively. They successfully reproduce 
the observed reorientation of P from P\\c to Pjja as 
a flop of the spin chirality from x||a to x||c. We 
determine the transition points and the spin struc- 
tures by calculating the T dependence of specific heat 
C{T) = j^d{n)/d{kBT) and spin chiralities XtiT) = 
7r(\Y.iiSi >^ S!.^f^)j\)/S'^ {l=a, b, c). Here the brackets 
denote thermal averages. Concerning the spin chiralities, 
the Xa{T) [xc{T)] has a large value, while other two com- 
ponents are nearly zero in the &c-plane [a6-plane] spiral 
or conical phases. Figure HJc) shows C{T) and X-y{T) at 
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FIG. 3: (color), (a) [(b)] Intensity map of internal magnetic 
field [Htn] for DyMnOa [TbMnOa] in plane of T and ex- 
ternal magnetic field H^^ [i^cx] calculated from experimen- 
tal magnetization data mb{T, H^^) [mc{T, H^^)], which re- 
produces the experimental T-He:^ diagram in Fig. me)[(f)]. 
(c)[(d)] Arrangement of the Dy [Tb] / moments in DyMnOa 
[TbMnOa] under H^^b [i/ex||c] where 0^60° [0~3O°] [|. 
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FIG. 4: (color online). Theoretical T-Hin phase diagrams of 
(a) TbMnOa and (b) DyMnOa for iJin||c, and spin structures 
in (c) PE(P=0) phase in TbMnOa and (d) FE(P||a) phase 
in DyMnOa. Dashed line (open circles) in (a) denote the 
crossover line (points) where C{T) shows a broad maximum. 



HP^=8 T for TbMnOs. The C{T) shows three peaks in 
accord with successive three phase transitions with lower- 
ing T. The assignments of spin structures are confirmed 
by calculating spin and spin-chirality correlations in the 
momentum space. 

By calculating the Hin dependence of the expectation 
value for each term in the Hamiltonian, we identify a 
mechanism of the chirality flop from xl|a to under 
Hcx\\b- The &c- plane spiral with x||a at iJex=0 is sta- 
bilized by the DM interaction associated with the DM 
vectors on the out-of-plane Mn-O-Mn bonds, i.e., 'Hf^M- 
The spins dominantly couple to the a-axis components 
of the vectors (i.e., components perpendicular to the be 
spiral plane) whose signs are the same within a plane 



but alternate along the c axis, while their magnitudes 
are all equal to ac [see Fig. [Hd)]. Without DM inter- 
action, angles between adjacent two spins along the c 
axis are uniformly 0c=7r because of the strong antifcr- 
romagnetic (ATM) coupling Jc- In the presence of DM 
interaction, the angles alternate between tt -|- A(j)c and 
TT - A(/)c with A0c>O [see the inset of Fig. [2jd)]. We 
can derive a gain of the DM energy due to this angle 
modulation as AE^^/N = -ac^^l cos (/)c| A0c- With- 
out Him the gain AE^^^ in the 6c-plane spiral domi- 
nates over the easy-(a6)-plane [or the hard-(c)-axis] spin 
anisotropy from 'Hf?^, which favors the a&-plane spiral 
with x||c. Note that the value of [cos 0c | is maximum 
(=1) at (/)c= TT, but decreases as 4>c decreases. This means 
that the application of Jfex||b suppresses this energy gain 
since it destroys the interplane AFM coupling and re- 
duces the angle (jjc from tt. The 6c-plane spiral becomes 
destabilized when the reduced energy gain AE^^ is de- 
feated by the easy-(a6)-plane anisotropy "H^g^, resulting in 
the spiral-plane (chirality) flop from be (xl|i) to a5 (x||c). 
Note that in i?Mn03, the ac-plane spiral or conical is 
unfavorable. This is because it can energetically benefit 
neither from nor from H^f^j , whereas the ab- and fee- 
plane spirals can take advantage of one of these two. We 
expect that the above mechanism is relevant also to the 
ffoxllb induced P fiop from P\\c to P\\a in LiCu202 
in terms of the role of Hex, which destabilizes the spin 
spiral with -P||c through destroying the AFM cou- 
pling along c. Note that the single-ion anisotropy H^i^ 
cannot work in this quantum S=l/2 spin system in con- 
trast to i?Mn03 with iS'=2 spins. We expect that the 
spin spiral with -P||a under H^x (possibly the a6-plane 
spiral) is stabilized by the other interaction, and the DM 
coupling with the c-axis components of DM vectors is a 
possible candidate. 

Now wc compare our results with experimental ones. 
Between Figs.[lje) and Fig.[2{b), there are a few discrep- 
ancies. First, threshold fields for the P reorientation are 
different; i.e., the calculated threshold value of H^^ for 
DyMnOs is approximately 18 T, whereas the experimen- 
tal value of Hex is 1-4 T. Second, the slope of the phase 
boundary is very steep in the theoretical T-Hm diagram 
of Fig. [2lb), while in the experimental T-H^x diagram of 
Fig. [He), it is rather gradual. These discrepancies are 
solved by considering the effective magnetic field Hfd 
generated by the rare-earth / moments, which acts on the 
Mn spins via the f-d coupling Jfd- Because of the AFM 
J/d, Hfd and H^x arc antiparallel, and the internal field 



Hi (7=a, 6, c) is given by Hi = H2x - H]^. Here H 



written using the /-electron magnetization m-y as a func- 
tion of T and H2x as Hj^{T,H2x) = zJfdrn^[T,H2x)- 
Here ^(=8) is the coordination number of R ions around 
the Mn ion. We assume J/d=0.45 T/hb for DyMnOa. 
Figure [3l[a) displays a color plot of the internal magnetic 
field Hi in the T-H^^ plane calculated using the experi- 
mental magnetization data. A solid line on which Hi is 
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equal to the calculated threshold value is drawn. This fig- 
ure coincides with the experimental diagram of DyMnOa 
in Fig. [TJc). A similar analysis for TbMnOa has also 
reproduced the experimental diagram (not shown). The 
roles of the f-d coupling in i?Mn03 at Hcx=0 have been 
studied by recent neutron-scattering experiments [l5l . [l6j . 
We find that the switching of P can be qualitatively un- 
derstood even without considering the f-d coupling, but 
it should be taken into account for quantitative discus- 
sion. 

Next we discuss the case of i?ox||c. The theoretical 
T-iJ^ phase diagrams of TbMnOs and DyMnOa are dis- 
played in Figs, mja) and (b). In Fig. IUJa), we find the 
transition to a coplanar spin state with P=0 for TbMnOa 
at i?^~3-5 T, which coincides with the experimental ob- 
servation of paraelectric phase under ffcx || c. For its mag- 
netic structure, see Fig. SJ^c). Again, there are a few 
discrepancies between the theoretical and experimental 
results [compare Figs. [Ijf) and Fig. |4|a)]. They are re- 
solved by considering the influence of Tb / moments. 
In Fig. [3Jb), we display the T and H^^ dependence of 
the internal field H^^ calculated from the experimental 
magnetization data. Here we assume J/d=0.65 T/fJ-B 
for TbMnOs. Solid lines on which H^^ is equal to the 
calculated threshold value are drawn. This figure coin- 
cides well with the experimental diagram of TbMnOa in 
Fig. [TJf). On the other hand, the transition to the ab- 
plane transverse conical state with P\\a [sec Fig. lUJd)] 
is found for DyMnOs in Fig. HJb) , which has not been 
observed in experiments up to H^y^=9 T. The required 
H^^ for this transition deviates from the calculated criti- 
cal value of iJfjj by the field Hfd from the Dy / moments 
antiparallel to iJex- Hopefully, the reorientation of P 
will be observed in DyMnOs under a higher H^^. 

The contrasting behaviors of P under ffoxljc between 
DyMnOa and TbMnOs can be attributed to the dif- 
ference in magnitude of the in-plane spin-exchange Jf,. 
TbMnOa has much smaller Jf,=0.64 meV than DyMnOa 
with Jh=0.99 meV. At ffox=0, the Mn spins form a spi- 
ral order to minimize the spin-exchange energy in both 
compounds. Once we apply i?cx||c, the ferromagnetic 
moment is induced along the c axis, and hence rotat- 
ing components of the spins become reduced. Then in 
TbMnOa with a small Jb, the spiral and conical spin or- 
ders no longer take advantage of the spin exchanges under 
ffcx II c, resulting in the first-order transition to the copla- 
nar state as shown in Fig. HJ^c). This state can benefit 
from all of the large a-axis components of the DM vectors 
on the out-of-plane bonds, which are perpendicular to the 
coplanar spin plane. The /fox-induced ferroelectric-to- 
paraelectric transition with sudden vanishing of P has 
also been observed in many other spin-spiral multifcr- 
roics, e.g., Ni3V2 08 and MnW04 We expect that 
the above mechanism is relevant also to them. 

In summary, we have theoretically studied the puz- 
zling T-ifcx phase diagrams of the spin-spiral multifcr- 



roic i?Mn03 (i?=Tb and Dy) and have revealed new 
mechanisms for the magnetic control of P by analyzing 
a microscopic spin model using the MC technique. We 
have shown that the applied HcxWQ {\\b in the present 
case) reduces the DM energy through modulating the 
interplane spin angles, and thereby controls a compe- 
tition between H^jyj and other interaction (H^j^ in the 
present case), which results in the spiral-plane or spin- 
chirality flop with reorientation of P. On the other hand, 
the applied Hcx^Q (||c in the present case) suppresses 
the spin-exchange energy through reducing the rotating 
components of spins, and thereby causes a competition 
between the spin exchanges "Hex and other interaction 
(^DM in the present case). As a result, the flrst-ordcr 
transition from spiral to coplanar spin phases occurs in 
TbMnOa with a rather small Jf, accompanied by the sud- 
den disappearance of P. We have discussed that the pro- 
posed mechanisms are also applicable to many other spin- 
spiral multiferroics. Additionally, we have found that the 
experimental results can be quantitatively reproduced by 
considering the effective field Hfd from the rare-earth / 
moments. 
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